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Abstract 

We consider the dynamics of a massless scalar field with time-dependent sources in the 
adiabatic limit. This is an example of an adiabatic problem without spectral gap. The main 
goal of our paper is to illustrate the difference between the error of the adiabatic approximation 
and the concept of non-adiabatic transitions for gapless systems. In our example the non- 
adiabatic transitions correspond to emission of free bosons, while the error of the adiabatic 
approximation is dominated by a velocity-dependent deformation of the ground state of the 
field. In order to capture these concepts precisely, we show how to construct super-adiabatic 
approximations for a gapless system. 

1 Introduction 

The adiabatic theorem of quantum mechanics is usually stated and proved for systems with a finite 
spectral gap. Mathematically it is known since at least 1998 O |3] that in a weaker sense the 
theorem remains valid even for systems without spectral gap as long as the spectral projections of 
the time-dependent Hamiltonian are sufficiently regular functions of time, see also [TZj . Its validity 
has been shown even for resonances [T] and for open systems [2]. In recent years the problem of 
adiabaticity for gapless systems attracted also interest in physics. For example, in [13] the authors 
consider adiabatic changes in the coupling of a scalar field. In a finite volume the ground state of 
the field is separated by a gap from the excited states. However, in the thermodynamic limit the 
gap closes and the spectrum of the field becomes absolutely continuous. The authors find that in 
the thermodynamic limit the error in the adiabatic approximation depends on the spatial dimension 
and other parameters controlling the spectral density of the field. They distinguish three regimes: 
one where the adiabatic theorem holds with the same error estimate as in the gapped case, one 
where it holds with a different power law and one where it fails completely. In all cases, however, 
they identify the error of the adiabatic approximation with the size of the non-adiabatic transitions. 
One main goal of our paper is to explain why this is, in general, not a valid identification and why 
the concept of super-adiabatic approximations is useful also in the gapless case. 

To this end we consider a specific gapless model in the adiabatic limit, namely a massless scalar 
field with time dependent sources. The goal is to exhibit a number of subtleties in adiabatic theory 
that are known in the gapped case - but often not expressed sufficiently clearly - also for the case 
of gapless systems. In particular we emphasize the difference between the error in the adiabatic 
approximation and what physically should be considered non-adiabatic transitions. To be more 
specific, let us briefly recall the situation in the presence of a spectral gap. Let H{t) be a time- 
dependent Hamiltonian and E{t) an eigenvalue with spectral projection P(t). Under appropriate 
regularity conditions on H{t) and the assumption that E{t) is separated by a gap from the rest of 
the spectrum of H{t), the adiabatic theorem states that the solution of the Schrodinger equation 

ief^^it) ^ H{t)il){t) with V(0)eranP(0) 



satisfies 

\\{l-P{t))m\\^0{e), (1) 



i.e. that any solution starting in the spectral subspace ranP(O) evolves into the subspace ranP(i) 
up to an error of order e. This error estimate, which contributes to the error of the adiabatic 
approximation, is optimal in the sense that generically the piece (1 — P(t))ip(t) of the solution in 
the orthogonal complement of ran P(t) is really of order e and not smaller for any finite time t where 
^H{t) ^ 0. However, if H{t) is constant outside of a bounded interval, say [0, T], and N + 1-times 
continuously differentiable, then clearly ip{t) £ ranP(i) for t < 0, but also 

P{t))tlj{t)\\=0{e^) for t>T. (2) 

Thus the non-adiabatic transitions into the orthogonal complement of ranP(t) after a compactly 
supported adiabatic change of the Hamiltonian are much smaller than 0{e). What happens is that 
during the adiabatic change of the Hamiltonian the solution does not exactly follow the spectral 
subspace ran P(t), which we will call the adiabatic subspace in the following, but slightly tilted 
e-dependent subspaces ranP^(i), so-called superadiabatic subspaces. With respect to the superadi- 
abatic subspaces one has 

P'{t))jlj{t)\\=0[e^) foralHsM. (3) 

In addition it holds that 

\\P^{t)-P[t)\\=0{e) 

which together with ([s]) implies the usual adiabatic estimate ([!]). Moreover, at times t when 
all derivatives of H{t) vanish, the superadiabatic projections agree with the adiabatic projection, 
P^{t) — P{t), and thus ^ implies also ([2]). In Figure [l] the situation is plotted schematically. 
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Figure 1: Schematic plot of the transition histories with respect to the adiabatic (full line) and the 
superadiabatic (dashed line) subspaces. The Hamiltonian varies only during the time interval [0, T] and 
outside of this interval the adiabatic subspace P{t), i.e. the instantaneous spectral subspace, agrees with the 
superadiabatic subspace P^{t). While the non-adiabatic transitions for t>T are of order e^, the norm of 
the piece of the solution that leaks into the orthogonal complement of P{t) for < i < T is typically much 
larger, namely of order e. The norm of the piece of the solution that leaks into the orthogonal complement 
of the superadiabatic subspace P'^{t) for < t < T remains small, i.e. of order e^, for all times. 



In this sense we say that the error of the adiabatic approximation is of order e, while the non- 
adiabatic transitions are (at most) of order e^. For a detailed exposition of adiabatic theory with 
spectral gap we refer to [TR' and references therein. Here let us mention only that the notion 
"superadiabatic" was termed by M. Berry in j4j- There he shows how for Hamiltonians depending 
analytically on time the transitions between superadiabatic subspaces are even exponentially small 
in e (as it is known also from the famous Landau-Zener model) and moreover, that the transitions 
as a function of time follow the universal shape of an error function, i.e. an integrated Gaussian. A 
rigorous account of [4 is given in [3 [Zl IB] ■ 
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In the present paper we show within a physically relevant but relatively simple example that in 
the gapless case the situation is to some extend similar. Before presenting our results in detail let 
us introduce the model and explain some of its important features. 

We consider a scalar massless field in three spatial dimensions. The momentum space for a 
single boson is and the state space of the field is the symmetric Fock space 

So '0 £ is a sequence ip = (i/'O: V'l: V'2, • • ■) with tpn{ki, . . . , fcn) a square integrable symmetric 
function of n variables in M^. For £ J- the inner product is 

The Hamiltonian of the free field is Hf dr(|/c|) and acts as 

{H{i>)niki, . . . , /c„) = J2i=i \ki\ipniki, . . . , fc„) . 

As a multiplication operator H{ is self-adjoint on its dense maximal domain D{Hf). Moreover, 
Hf has a unique ground state given by the Fock vacuum :— (1, 0, 0, . . .). Now we add moving 
charges as sources to the field. For notational simplicity we assume that all sources have the same 
normalized form factor : M'^ — > [0, oo) with 

/r3 dx (p{x) = 1 . 

Later on we will make additional assumptions on the "charge distribution" ip. The sources are 
located at positions Xj £ M.^ with total charge ej and form factor ip. The linear coupling to the 
field is given by the operator 



|fe| 

where we abbreviate x = (xi, . . . , X]^). Here if denotes the Fourier transform of ip and $ is the field 
operator 

$(/) :=^(at(/) + a(/)), 
where a and are the standard bosonic annihilation and creation operators on J-', see Section 



3.1 



for more details. Under appropriate conditions on ip the total Hamiltonian H{x) = Hf -\- Hi{x) 
is self-adjoint on D{H) :— D{Hf) and bounded from below. If the total charge of the sources 
e — is zero, H{x) has a unique ground state ^ JIq- The ground state U{x) contains 

so-called "virtual bosons" that provide a "dressing" of the sources. In contrast to "free bosons", 
the virtual bosons do not propagate. These heuristic notions will be made more precise later on. If 
the total charge is different from zero there is still a good notion of approximate ground state, at 
least for our purposes. 

We will be interested in the dynamics of the scalar field when the sources change their locations 
with time. Let : M — )• M.^, 1 1— Xj{t) be smooth, then 

H{t) := H{x{t)) = H, + Ef=i Hi^A^,{t)) 

describes the evolution of the scalar field in the presence of sources moving along the prescribed 
trajectories Xj{t). We consider the solutions of the Schrodinger equation 

ief^m = H{t)m 

for the field in the adiabatic limit e <C 1, which corresponds to slowly moving sources. Assume for 
a moment that the field is restricted to a finite box. Then the Hamiltonian H{t) has purely discrete 
spectrum and the ground state is separated from the first excited state by a gap. The adiabatic 
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theorem implies that if the field starts in the ground state il{x{0)) of H{0), it remains up to a 
phase e-close to the ground state V,{x{t)) of H{t) at later times. The solution is thus e-close to a 
state where a static cloud of virtual bosons follows the slowly moving sources and no radiation is 
emitted. The error of the adiabatic approximation is of order e. If we pass to the superadiabatic 
approximation, we realize that the solution is even e^-close to a state W^{t) where an £-dependent 
cloud of virtual bosons follows the slowly moving sources and no radiation is emitted. The state 
ri^(t) is still e-close to the static ground state however, the dressing by virtual bosons 

depends now also on the velocities Xj (t) and higher derivatives of Xj (t) . If at some time the sources 
come to rest then from that time on 17^ (i) = ^l{x{t)), hence the field is e^-close to the static ground 
state again. In particular, the probability for emitting a free boson is at most of order e'^^ . 
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Figure 2: The situation with spectral gap; For t < the sources are at rest and the adiabatic approximation 
agrees with the true solution. While the sources are moving, the dressing is deformed, which leads to an 
error of order e in the adiabatic approximation. The superadiabatic approximation takes into account the 
deformation and the error is only of order . When the sources are again at rest for t > T, the dressing 
is again the static one and the error of the adiabatic and the superadiabatic approximation is of order e^. 
This error is due to non-adiabatic transitions that correspond to the emission of free bosons. 

Summarizing the above, we have that as long as there is a gap and as long as the prescribed 
trajectories Xj{t) are C^"'"^(M, E'^), non-adiabatic transitions corresponding to the emission of free 
bosons are at most of order . Nonetheless, the error of the adiabatic approximation is of order 
e for all times where Xj{t) ^ because of the deformation of the cloud of virtual bosons, see also 
Figure [2] This exemplifies in a simple physical example why it is important to distinguish carefully 
between the error of the adiabatic approximation and the non-adiabatic transitions. It also is 
an example where superadiabatic subspaces have a clear physical meaning as velocity-dependent 
dressing of the sources. This is of course all well known and the content of this paper is to show 
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that this picture survives also in the case without spectral gap, i.e. in the thermodynamic limit 
where the box is replaced by K^. 

In the remainder of the introduction we provide informal statements of our main results, where 
for simplicity we restrict to the case of a neutral system. Then there is a unitary transformation 
V{x) e ^{J-) such that V{x)H{x)V{x)* = Hf + E{x) and thus H{x) has a unique ground state 
U,{x) = V*{x)Q,Q with eigenvalue E{x). The general and rigorous statements are explained in 
Section 2. Our first result is an adiabatic theorem without spectral gap, cf. Theorem [ij 

Adiabatic Theorem: The solution of 

ief^m = H{x{t))4,{t) (4) 
with V'(O) = r2(a;(0)) satisfies for any < G M 

||^(t)_e-i/od«^(-(-))l7(2;(t))|| ^©(eT^TTi)) . (5) 

Thus, although the spectrum of H{x) is the whole half line [E{x), oo) and the eigenvalue E{x) 
is not separated by a gap from the continuous spectrum, the usual adiabatic approximation for 
the evolution of eigenstates holds with almost the same error estimate as in the gapped case. We 
have even more: if one adds free bosons then their dynamics is adiabatic too in the following sense. 
To each configuration x of the sources the annihilation and creation operators of free bosons are 
V{x)*aV{x) and bl := V{x)*aW{x) and thus 

ei^(-)*6f(/)e-i^(-)* = 6f(/(t)), 

where f{k,t) — e^el'=l*/(fc) is just the free time evolution of a boson and fof" stands either for or 
Now we define the adiabatic approximation as 

m 

■.= e-ifo<^'^E^'^i'^))Y[bl^^^{fi{t))n{x{t)). 
1=1 

We will show that it approximates the solution of ^ with initial datum 

1=1 

in the same sense as in ([s]), 

\m)-i^.d{t)\\^0{eVHm)- (6) 

If we denote by Qm{x) := ^^(2^)* Qm the projection on the sector of Fock space containing 
exactly m e No free bosons, where Qm is the projection on the m-particle sector of Fock space, the 
above result implies that for ip{Q) E Q,n(x{Q))J^ 

- QmixmmW = 0{e^Hl/e)) , 

i.e. that the probability for emitting a free boson is at most of order ln(l/e). Note that this means 
that not only the spectral subspace QQ{x(t))J- is adiabatically invariant, but also the subspaces 
Qm{x{t))J- , which are not spectral subspaces of the instantaneous Hamiltonian H{x{t)). 

In order to understand what part of (1 — Qm{x{t))) i){t) really corresponds to emission of free 
bosons and what part is merely a velocity-dependent deformation of the dressing, we need to 
introduce the superadiabatic picture. Now not only the dressing of the sources depends on their 
velocities, but also the annihilation and creation operators of free bosons and the associated m- 
particle sectors of Fock space. In Lemma |9] we construct a velocity-dependent dressing operator 
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V^{t) such that the corresponding objects n%t) := V*{t)no, := V^* {t)a*V§{t) and (5^(t) := 
V^*{t)QmV^{t) allow for an improved adiabatic approximation. Let 

m 

1=1 

be the superadiabatic approximation to the solution ■i/;(t) of Q with initial datum 

m 

1=1 

then ^ can be improved in the following sense: the part of the solution tp{t) that stays in the 
TO-particle sector is given by V'su(i) with a much smaller error, 

\\Q'„Mm-i'snm = o{sHnii/£)). 

However, more importantly, the superadiabatic representation can be used to compute the non- 
adiabatic transitions corresponding to the emission of free bosons by first order perturbation theory. 
In Theorem [2] we show that the "non- adiabatic" wave function 

^naW -ie^*dse-j^= (e-^l^'lVW) Auis) 

satisfies 

||(l-Qf„(t))^(0-V'na(<)|| =0(e'ln(l/£)) . 

Here '■— ^ (&s^(/) + ^s(/)) ^'^'^ 9'^{s) — 9'^ x{s)) is a couphng function depending on 

the acceleration of the sources. In summary 

\\m - (V-suW + ^n.m\ = 0{e^ ln(l/e)) 

and thus '4'suit)+ tpna.it) yields a good approximation of the true solution ip(t) with a clear separation 
of a "superadiabatic" piece ipsu{t) that improves the adiabatic approximation and a "non-adiabatic" 
piece ipnait) that contains real non-adiabatic transitions. Note that a posteriori one can compute 
the leading order of ipnuit) replacing all superadiabatic objects by their leading order adiabatic 
counterparts, 

^„a(t) = -ie£dse-?^*'^^'^(^(^~» {e-'^'^'g'is)) + 0{e^Hl/e)) . 

Finally one can ask for the probability of emitting a free boson, i.e. for computing HV'nalOlP- 
However, due to the infrared problem it turns out that this has no nice asymptotics for e ^ and 
we can only show HV'naC^)!!^ = 0{\t\'^e^\n{l/e)). But a physically more interesting question is the 
energy radiated through the emission of free bosons. In Theorem [3] we show that when starting in 
the dressed vacuum, i.e. ip{0) — i^'^(O), then 

E,,d{t) (V^„a(t),i?f(t)V'na(t))^ = 2^ ^ ds \d{s)\^ + o{e^) , 

where 

is the second derivative of the dipole moment of the sources and H{{t) is the field energy of the free 
photons. Taking a derivative gives the Larmor formula for the radiated power of slowly moving 
charges, 
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To come back to our main message once more, note that computing naively the energy in the 
piece of the solution that constitutes the error of the adiabatic approximation ^{t) := ijj{t) — V'ad(0 
would have led to 

where k :— k/\k\. As our analysis showed, this velocity-dependent term is the energy in the 
deformation of the dressing and being of order e^, it is much bigger than the energy emitted by 
radiation of free bosons, which is of order e'^. In order to obtain the correct picture, the use of 
superadiabatic approximations seems inevitable. 

Finally let us mention that a closely related problem was considered in [TB] , see also [13] • There 
semiclassical limit of non-relativistic Schrodinger particles coupled to a massless scalar field was 
considered. In a sense the system considered in the present paper can be considered a toy model 
version of the one in |16| . However, as a consequence of the simplicity of the model the results 
obtained here are much stronger and at the same time the proofs are more transparent. For 

Acknowledgements. We thank Rainer Nagel for pointing out to us reference [llj . 



2 Main results 

We always assume the following properties for the parameters of the model: 

Assumptions. The charge distribution ip G S'(M'^) is a spherically symmetric Schwartz- function 
with ^(0) = (27r)~3/2 ^. g (^4^]g^ ]g3) /or j = 1, . . . , N. 

As explained in the introduction, the central object that we construct is the transformation to the 
superadiabatic representation, a unitary dressing transformation 

V'{t) -.T . 

In the superadiabatic representation the Schrodinger equation reads 

V-(t) (ieA _ ii{x{t))) V-(tr V%t)m =: (ie^ - . 

In this representation the Fock vacuum f2o corresponds to the dressed ground state and the m- 
particle sectors of Fock space correspond to states with m free bosons. While in the introduction 
we formulated the statements in the original representation, it is more convenient to formulate 
and prove everything in the new representation after performing the dressing transformation. The 
translation back to the original representation is straightforward using the definitions of fof"^ of the 
introduction. Note, however, that if the total charge is different from zero, the transformation ^^(i) 
has no asymptotic expansion in powers of e, not even the limit lim£_j.o V^it) exists. 

Let Q<M denote the projection on the sectors of Fock space with at most M photons. 

Theorem 1 (Adiabatic approximation). The dressed Hamiltonian H^^^^^ generates a unitary 
propagator ?7drcss(i) C^drcssC^j^o = 0). For any T £ M. there are constants C < oo and Eq > 
such that for all Af G N, < e < and \t\ < T 

in the norm of C{J-) and in the norm of C{D{Hf)). Here 

E{t) :^-|43dfc (7) 
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with 

-o(-W,fc):=E|pe--^(*). (8) 

Thus at leading order the time evolution is just the free evolution of the bosons and a time- 
dependent dynamical phase. In particular, the m-boson sectors QmJ' are adiabatic invariants of 
the dynamics. An initial state in Qm^ of the form 

evolves into 

Applying V'^{t) and using 

\\V%t)-Vixm^Ois) 
in the case of a neutral system yields 

Theorem 2 (Superadiabatic approximation). For any T e M there are constants C < oo and 
eo > such that for all Af G N, < e < eo and \t\ < T 



(u,,,Ut) - e-i(^'*+/o ^^(^)'i^) (^1 + i£^*d.e*^'f /...ad(s)e-'^'f)) 



Q 



<M 



< 



< Ce^ ln(£-i) VmTi 

in L{F) and C{D{H[)) with 

E'it) = E{t) + f E / dfc ?Z^lM(^e-"=-(-W-^(*)) {n ■ i,{t)) ■ x^) 
4 j^i-^R^ 1^1 

and 

N 



it) :=cI>(l[,.,^)(|fc|)^^e->'=-^(*)(.,x,(t))) 



Note the infrared regularization by the characteristic function Ij^s oo)(|fc|) cutting off boson 
momenta smaller than in this definition, which can be omitted for neutral systems. 

With the help of this theorem we can define a decomposition of the wave function into a supera- 
diabatic part and the non-adiabatic part 0„a with different photon number. For a wavefunction 
(/'o G QmJ' we define 

and 



Our final result concerns the amount of energy lost by the system due to radiation. For the 
stationary problem with Hamiltonian H{x) independent of time, the natural definition for the 
energy carried by free photons in the state ^ e D{H[) is 

-^rad.stat 

where E(x) is the ground state energy. For time-dependent Hamiltonians the definition of energy is 
somewhat subtle, more so, because we look at the energy of a subsystem, the free photons. However, 
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in the superadiabatic representation there is again a natural definition for the energy given by the 
free photons, namely 



Er^Mit)) := {V'{t)tP,HfV'{t)ij). (9) 
We will show that at any time t when x{t) — x{t) — 0, it holds that 

^rad(l^(t)) = ^rad.statlV'Ct)) + 0{e^) . (10) 

Thus in a situation where the sources move only during a finite time interval, before and after the 
change the definition of energy agrees with the static one. Note that the error of order comes from 
the infrared regularization in the transformation V'^(t) and is identically zero for neutral systems. 
The equality ( 10 ) is also determining the exponent 8 in the infrared regularization. A smaller 
exponent there would lead to a larger error in (10). 

If we assume that the initial state is the dressed vacuum, then the energy emitted as a result of 
the motion of the sources is just the energy of the free photons and thus ii'rad(^) — (0na(i), ^^f0na(O) ' 
which is computed explicitly in the following theorem. 

Theorem 3 (Radiation). Let 0(0) = ^l, then uniformly on hounded intervals in time it holds that 

/o ^0 



E^Ut) = y E ^.^^/ '^^ll^/ / d.d.'e'l'=l^e-"=-(-^(^)--(^')) (11) 
X K ■ ij{s) K ■ i,(s') + O (£''(ln(£-i))^) 
ds|Q!'(s)p + o(e3), (12) 



£3 



2-67r 



where 



is the second derivative of the dipole moment. 



3 Mathematical details 
3.1 Operators on Fock space 

In this section wc recall the definitions of the operators on Fock space and some of their important 
properties. Proofs of all claims can be found in [14 , Section X.7. 

We call J'fin the subspace of the Fock space for which ipm = for all but finitely many m. 
The second quantization of a self-adjoint multiplication operator tu with domain D{w) e L^(M"^) is 
defined for i(j e J^fin with components ■0™ S ^T^i-^i"^) ^ 

(dr(w)?/')m(fcl, • • ■ , hn) = J2T=1 ^{kj)i^m{kl, . . . , fc„) 

and is essentially self-adjoint. In particular, the free field Hamiltonian iJf :— dr(|fc|) is self-adjoint 
on its maximal domain. The annihilation operator and the creation operator on ■!/; G J-fin are defined 

by 

(a(/)?/')m(fcl, ■ • ■ , fcm) V"l + 1 /jj3 dfc/(fc)V'„j+i(fc,fci, . . . , k,n) , 

{a\f)'^)m{ki, . . . , fc™) := ^= Ylj-i f(kj)ipm-i{ki, ...,kj,..., fc„) . 

Here kj means that kj is omitted. They fulfill the canonical commutation relations 
K/),a^(3)] = (/,5)L^(M3), [a(/),a(5)]=0, [a\f),a\g)]^0. 
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The creation and annihilation operator define the Segal field operator 

It is essentially self-adjoint on J-fin- The canonical commutation relations imply 

[$(/), a>(.g)] =ilm(/,5)i2(K3) and [dFH, i<i>(/)] = <i>(ic^/) . (13) 

3.2 The Hamiltonian 

As described in the introduction we consider the Hamiltonian 

H{t) := H{x{t)) - Hi + $(«o(a:(t), fc)) Hf + $o(0 

with vt^{x{t), k) defined in ([S]). As to be stated more precisely in Lemma[6]and after, this Hamiltonian 
has a ground state only for a neutral system, i.e. if ~ ^^ otherwise the bottom of the 

spectrum is not an eigenvalue. To deal with this fact we introduce an infrared cutoff < cr < 1 in 
the interaction and put ^Pa-ik) = l[^^oo)i\k\)V'{k) and v„{x,k) — 1^^,00) (|fc|)^^o(2;, fc)- The resulting 
truncated Hamiltonian {t) — H{ + (t) is a good approximation to the original Hamiltonian as 
will be proven in Proposition 11 We emphasize that H"'{t) is used only as a tool in the proofs and 
the final results of Section 2 hold for the Hamiltonian H{t) = H'^^'-\t) without infrared cutoff. 

Lemma 4. The operators H'^{t) are self-adjoint on D{H{) for all t and a € [0,1]. The graph 
norms of {t) are all equivalent to the one defined by Hf uniformly in a. 

Proof. There exists a standard estimate for $ (cf. e.g. [5 Proposition 1.3.8). Let e D{Hf) and 
6 > 0, then 

WHfHWr < Mi^f^ii^ + (il|//^lH.(K3) +2ii/iii.(K3)) ii^ii^ . 

Thus $(/) is infinitesimally bounded with respect to H{ if 



J3) + ||//v4fc|lli2(R-') < 00. (14) 



Equation (14) holds uniformly for / = Vo-{x{t),k) and a G [0,1]. Hence H"'(t) is self-adjoint on 
D{H() by the Kato-Rellich Theorem (cf. [14j, Theorem X.12) for all a S [0, 1] and the graph norms 
are all equivalent. □ 



Proposition 5. Let t 1— > x{t) e C"(M,M'' ), n £ N, and equip D{H{) with the graph norm. Let 
(7 G [0, 1], then 

2. H'^ generates a strongly continuous unitary evolution family denoted by UH'{t,ti:i)- Moreover, 
i7//<7(i,to) is a bounded operator on D{H{) with 

sup ||?7fl--(s,to)|l£(c(Hf)) < 00 (15) 
se[to,t] 

for any t G M and uniformly in a. 



Proof. According to (|14|) we need to show that 

< 00 . 



I dt" 



vo{x{t),k) + ^vo{x{t),k)/^\k\ 



which obviously holds under our assumptions. That H and H" generate unitary evolution fami- 
lies follows from the first statement and the general result about contraction semigroups (e.g. [14j . 
Theorem X.70). The norm bound Uwit^to) in C{D(Hf)) is less known and was proved in Theo- 
rem 5.1.(e) of dl]. □ 
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Lemma 6. Let K = H[ + $(z(fc)) with z{k) e L^(R^) such that 

" 2 Jr3 |fe| °o ■ 

Then the spectrum of K is given by [a, oo) and the infimum of the spectrum a is an eigenvalue if 
and only if 

In this case the unitary operator V = e^'*''^*^'^)/''^!) is well-defined and K = V*H{V + a. 

Proof. The first part is Proposition 3.10 of [TU] and the second part is Proposition 3.13 of the same 
paper. □ 

Corollary 7. The spectrum of Hit) is cr{H{t)) = a^c{H{t)) = [E{t), oo) with E{t) as in 

Corollary 8. For a > 0, the infrared regularized Hamiltonian H"{t) can be written as 

H^t)^V:it)HtV^{t)+E,{t), 

where 

Its only eigenvalue with eigenvector Qa{t) ■= V*{t)flQ is 

E,{t) = E{t) + 0{n) . (16) 
Proof. The last statement is an immediate consequence of Lemma [6] and the fact that 

\v,{x{t),k)\^ f \vo{x{t),k)\' 



dfc ' , ' " = / dfc ' ' " +0{a). □ 

\k\ Jr3 |fc| 

3.3 Superadiabatic perturbation and the dressing operator 

In this section we introduce the perturbed projections and the transformation which "diagonalizes" 
the infrared regularized Hamiltonian up to 0{s'^). As explained in the introduction the idea is to 
modify the adiabatic projections PS^it) := V* {t)QmVo-(t) on the m-free-photons subspaces in such 
a way, that the modified superadiabatic projections P^'^(t) are higher order adiabatic invariants, 
i.e. that 

[iei-^~H(t),Pr^\t)] =0{s') 
holds. This construction is by now standard, see |12| . and yields 

= V:{t)QMt)+ieV:{t)[Q,^,<f2{t)]V„{t)+0{e^) , 

where we introduce the shorthands 



^,it) := -i[dr(|fcri),<i>i] = $(i^ e'^-^^^'\^,x,it))) ^z,it)) . 



\k\ 

Here k is the unit vector in the direction of k. The canonical commutation relations link $i and 
$2 in the following way 

i[dr(|fc|),$2] = $1. 
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Now the corresponding dressing operator V^{t) that maps P™'^(i) to Qm then needs to have the 
expansion 

V^{t) = (1 + ie$2W + 0{e^))V^{t) and V^{ty = y;(<)(l - i£$2 W + 0(£')) . 
This suggests to define 

Lemma 9, For a G (0,1] and e £ [0,1] the operator V^{t) is unitary, belongs to C{D(Hf)) and 
SQ/tzsJi/es 

wv^mc {D{Hf)) + \\Kity\\c{D{Hi)) <C <00 

uniformly in a and e. 

The map M — > C{D(Hf),T), 1 1— > V^{t) is differ entiahle and 

f^V^ = (-i$i+ie$2+£[$2,$l]-4[*2,*2]) Kx'- (17) 



Proof. With 



and 



[iJf,K(t)] = V„{t){V:{t)HiV,{t)~Hi) = y,(t)($o(t)-Sa(i)) 



f 

[-i^2,hi$2,hi^2,---gfjjjj _ 

e-*^ (e$i-i4[$2,<fi] 



the first statement foUows, since $o, 'i'l and [$2, 'I'l] = ilm(z2, zi) are bounded independently of a 
as operators from D{Hf) to T . 

Under our hypotheses on (p the operators $(iwcr(a;(i), fc)/|fc|) and $ = $1 commute as a conse- 
quence of (13l and {v„{x{t),k) l\k\,Va{x{t^,k)n ■ Xj{€)) =0. Hence 

= ^e-'*('-'(-(*)-'=)/l'=l) = -ici>(i)T4(t) = -i<J>i(t)F,(t). 
For the other exponential we find 

-le-*^ ^ [at,e'^*^] = (5t-e'^*=ate-'^*^)e'^*^ = (i£$2 - 4[*2, $2]) e'^*^ . 

Finally 



implies (17). □ 
3.4 The dressed Hamiltonian 

With the help of the dressing transformation V/(i), we define the dressed Hamiltonians 
ieat-i?drc.sW V^{t){iedt~li[t))V^{t)\ 

As a consequence of Lemma |9] the dressed Hamiltonians are self-adjoint on D{H{) since V^{t) is a 
bijection on D{Hf). The corresponding evolution families are related through 

Udr-Ut) = V^{t) Unm^iOy and [/^^.^(i) = V;{t) Uh^ m^iO)* . 



12 



Lemma 10. The dressed Hamiltonian has the form 
with 

El{t) := E,{t) - 4lni(z2W, ^iW) + ^Im(z2(t), W) ■ 
Proof. For H'' {t) we find 



ye jja ye* 



1=0 



[i<&2,[i'i>2,[i'i>2,---gf]]]] , ^ 

7. h -C/n- 



= iJf + -i^ [$2,4-1] = Hi + E„ - e<i>i + ^lm{z2, zi) . 



and with (17) 



FJieaty^ = ie9t-i£(3^^')^t^" = ieS^ - £$i + e^^^ - ie2[$2, $i] + i^[$2, $2 



Hence 



Hdressit) = Hi + E,{t)-e^<i>2{t)-'-^lm{z2{t),zi{t)) + '-^lm{z2it),Z2{t)). □ 
Proposition 11. For any a,e E (0, 1] and M e N we have 

I (c/drcssW - f/d"css(i))Q<M|| = O (n/M+1 ^) 

uniformly on bounded intervals in time in C{T) and C{D{Hf)) . 
Proof. By definition we have 

Let 

His) - H^s) = <f{^o,.m)vo{x{s),k)) =: <l>(/<.(s)) . 
Then bosons with momenta smaller than a evolve freely under the regularized evolution, i.e. 

{7^„(t)<i>(e'l'=li/<,(s))C/H.W = $(/<,(s)) 

holds for all t e M. To see this note that the equality holds for t = and the derivative of the left 
hand side vanishes since 

[i^^<f(e'l'=l^/<.(s))] = [i/f,<i>(c'l'^-lt/<,(s))] = -i£Aci,(ei|'=l!/<,(s)). 

Moreover, ^{f<a{s)) commutes with V^{0)* . By the Duhamel formula we thus find 

WiUHit) -UH^{t))V^{Or Q<m\\ < ^ fds \\{H{s) - H''{s))UH^{s)V^{0rQ<M\\ 



ds \mf<,{s))UH^{s)V^{0)*Q<M\\ 



ds 



C/ff.(s)$(-e'l''^lf/<„(s))F;(0)*Q<M 

- f ds $(-e'l'=lf/<,(s))Q<M 
£ Jo " " 



The following lemma, which will be applied several times in the following, together with the fact 
that 

\\f<.mL^^ + \\\k\f<As)\\L-^^Oia) 
uniformly for s € [0, t] allows us to conclude. □ 
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Lemma 12. There is a constant C < oo such that for any f E L'^(R^) and AI E N 

mf)Q<M\\c(^) < 25^/mTT||/IL.(r3) 



ll*(/)Q<M||£(D(ff,)) <CViWTI(||/|L.(R3) + |l|/c|/||^.(K3)^ 

Proof. According to [H], Theorem X.41, we have for / e L2(M^) that 

ll$(/)Q<Af^||^ < 25VmTi ||/IL.(R3) . 

The second claim follows from the observation that 

HMf)Q<M = Hf)Q<MHf - i$(i|fc|/)Q<Af 
together with the first estimate. □ 

3.5 Effective dynamics 

In this section we first show that the statements of Theorems [T] and [2] hold for the infrared regularized 



evolution C/dressC^) with an error depending only logarithmically on a. Then we use Proposition 11 
and an appropriate choice for a — cr{e) to show the statements also for the full evolution [/dress 
Clearly the statement of Theorem [l] is a consequence of Theorem [2j but since the proof necessarily 
proceeds in the same two steps, we separated the statements. 

Lemma 13. For any < a < ^ and M eN it holds that ^2Q<m G C{D{H{)) n C{J^) with 



||*2Q<A/|| = 0(VMTTln(a-i)) 
uniformly on bounded time intervals in C{J-) and C{D{Hf)) . 
Proof. Since 

||z2|li.(K3) = / dfc|f;'^e-*-^ (^(«:,£,) + (^,i,)2) ' = 0(ln(a"i)), (18) 

the claims follow from Lemma [T2l □ 
Proposition 14. For cr, e e (0, |] and M e N we have that 

uniformly on hounded time intervals in C{J-) and C{D(Hf)) . 

Proof. Note that (15) together with Lemma [9] impHes C/dross(^) ^ uniformly for s e [0,i]. 

With the shorthand Uo{t) := e--AHi^+!Q El(s)ds) g^^j ^-^^ 

previous lemma, we thus find 

C^dress(i)Q<M = Un{t)Q<M + -US,,ss{t) ^ ds U^.^Jis) $2(i) Uo{t) Q <M 

£ Jo 

= Uo{t)Q<M + i£ C/d"css(t) / ds C/d",.ess*(s) ^2{t) Q<M Uo{t) 

□ 







Uoit)Q<M+o(^s\n{a-')^ 
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For the next proposition we abbreviate 



with 



rad 



Proposition 15. For a,e E (0, ^] and M G N we have that 

I (c^drcssW - C^cffW)Q<M|| = VmTT In {a'')) 

uniformly on bounded time intervals in C{J-) and C{D{Hf)) . 

Proof. Using again the Duhamel argument together with Proposition [14] we find 

t/drc.s(^)Q<M = U,{t)Q<M + >d"css(0 ds C/d",_*(s) ^2{t) Uo{t) Q<M 



(19) 



(7o(t)Q<Af + ieC/drcss(i) / dsC/d",,,,*(s)Q<Af+i$2(i)0<MC/oW 
Jo 

(7o(t)Q<M + i£C/oW / dsU;{s)^2{t)Q<MUoit) + o(^e^Ha-') 



in the norm of £(I?(iJf )). Inspecting (18 1 shows that 

N 



$2 = /^s;d + * 



Hence ( 19 ) follows, once we show that the second term does not contribute to the leading order 
transitions. This follows from an integration by parts, 

^ dse'"'i^{fAs))e-'"'iQ<M = ^ ds$ (e-'l^^'lf Q<m 



\e i As $ 



|fc| 







^ Q<M 






Jo 



(^e-'I'^lf (£/.(.))) Q<M 

□ 



Now Theorem 1 and 2 follow from Proposition 15 and Proposition 11 To see this put cr(e) := e 
and define 

and U',^{t) -.^ U:^'\t) . 

Then 



(C/dress(i)-C/effW)Q<M 



< 



< 



C/dressW - C/dll(0)Q<Af|| + || (c/dHW " C^cffW)Q<M 



= ©(e^ Vm + 1) + ©(e^ In(e-i) VM + 1) , 
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3.6 Non-adiabatic transitions 



It is now evident how to define the part of the wave function that corresponds to emitted photons. 
As an immediate consequence of Theorem [2] we obtain the following corollary. 

Corollary 16. Let (/)(0) € QmD{Hf) and a{e) = then 



4>sn{t) := 
have the expansions 



,C/drcs.(i, 0)0(0) and 0„aW QmC^dross(t, 0)0(0) 

0,„(t)=e-i(^'-*+^*^^(^)^^V(O)+i?i 



0„a(t) = iee-^(^'*+^'^(^)'i^) / dse'"^ih,,d{s)e-'"^mO) + R2 



and 



with -0(£2ln(e-i)) and \\R2\\d(h,) = 0(e' Me-^)) . 

The probability for emitting a photon mitil time t is thus given by ||(/>na(i)|P- But this has no 
simple asymptotics for e — >■ because of the infrared problem. However, for the radiated energy 
there is a simple asymptotic expression. 

Recall the definition of the energy of free photons in ([9]). Note first that at any time t when 
x{i) = x(t) = 0, it holds that Vit) = V^(_^)(t) and E%t) = E{t) and thus with Corollary [s] 



In the third equality we used equation (14) to bound \\H — -ff'^||£(£)(/ff) jf) by Ca^/^ and (16). This 
shows that whenever the energy of free photons is unambiguously defined, our definition using the 
superadiabatic approximation agrees with it. 

We will now prove Theorem [3] by plugging the superadiabatic approximation into the definition 
([9]). Note that this computation of the energy is the reason that we insisted on all our estimates 
being vahd also in £{D{H{)). 



Proof of Theorem^ For (11 ) observe that with the splitting given in Corollary 16 we have that 

E,,dWt)) = {(t>{t),Hi4>{t)) = {(p,^{t) + 4>,Ut),Hi{(t)Ut) + 0na(t))) 

= (0na(O,^f'^na(t)) = ^?ad(i) + 2Re(i?2 , iff </>„a (0) + (^2,i?ti?2) , 

where £'rad(^) J^^^ stands for the explicit expression claimed in ( 11 ). For the last term {R2, HfR2) - 
0(e'' ln(e~^)^) follows from ||^2||£i(jjj) = C(e^ln(e~^)) and for the mixed term 

|(i?2,i?f0„a(O)l < \\R2\\^\\Hf<l>n.{t)y = o(^e\He-^)^)) 



from integrating Hf(j)^i^(t) by parts as in the proof of Proposition 11 



In order to get also ( 12 ), we first transform the integral (111 with r' 
a(r) := min( — , -{t — r)). The resulting integral is 



and define 



N 



Co Ci 



dr \(pair)\'' 



a(r) 



'0 J -a{T) 

Kr, Ujiir, £t', k) := k ■ Xj{T + ^) k ■ Xi{T 



dr dr' e'''(^'-''^'(^'^^'^'^»aj,(T, er', k) , 

) and f3ji{T,eT' , k) := 



where we abbreviated k = nr, aji{T, er' , k) := k ■ Xj{T + ^) n ■ Xi{T — c^iii^ 
K ■ {xj(T + ^) — Xiir — ^)). Note that aji and are uniformly bounded on the domain 
of integration. Thus we can replace |i^ct('')P by |(^(r)p with an error of order e^tTjip, which is 
negligible. 
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Observing that after summation the integrand is symmetric with respect to change of sign in 
r and r', i.e. for fj^{r,T') |^(r)|2e''^(^'-ft-'(^'^^''«))aji(r, er', k) it holds that Y.^,j f]^ir,T') = 

Tlii j fjii^'^'i ^■'''): can extend the r-integration to aU of M at a cost of a factor i. Modulo higher 
order terms this leads to 

N , ft /.a(T) foo 

\2Ar(T'-Pji(T,eT',ii)) 



I{t) = — 2, Sj^i / / / dr dr' aji(r, er', k) / dr 



|^(r)|^e" 



y] Cje, dK / / drdr'ajj(r,£T',K)(^|(^p)(r' - /3ji(r,£r',K)) , 

; JS^ Jo J-a(T) 



4 

where ^ denotes the Fourier transformation. Since with (p also is a Schwartz function and 

since aji and are uniformly bounded on the domain of integration, we can shrink a(T) to tq with 
an error that is asymptotically smaller then any inverse power of tq in the region where tq < a(r). 
The same replacement in the region where tq > a{T) leads to an error of order ©(e^ro). Taking e.g. 
To ~ E~^l'^ leads to negligible errors. Now we can Taylor expand the integrand as 

gir(r'-ft-.(r,5r'))^^.^(^^^^') = e''"^"'-'^^' ("'°)+°(^""» (a^, (t, 0) + O (eTo ) ) 

= e''-(^'-^-(^-"))a,,(T, 0) + + r)eTo)) 
and the remaining leading order contribution is 



I[t) — — y CjBi / dn dr K ■ Xj{T) K ■ Xi{T) / dr' / dr|<^( 

y CjCi / dn dr K ■ Xj{T) K ■ Xi{T) df dr\tp{r)\'^e 

, , n JS^ Jo J-fn J-oo 



7-)|2gi'"r'g-i'-K-(a:j(r)-a;i(T)) 



with f — t' — K ■ (xj^r) — Xi{T)). With the same argument as above we can replace Tq by oo while 
making an negligible error and end up with 

3 ^ r rt 



I{t) = — - — J" ejd dn dTK■XJ{T)K■it^{T) df(^|^n(f) 

N f* f 

'(0)P / ^j^i dK K ■ Xj{T) K ■ Xi{T) 

"'0 Js^ 



4 

27r£3 1 /"* 



4 (27r)3 3 



^ e^ei / dr Xj(r) • aji(T) . □ 
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